Quantizing magnetic field and quark-hadron phase transition in a neutron star 
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We investigate the influence of a strong magnetic field on various properties of neutron stars 
with quark-hadron phase transition. The one-gluon exchange contribution in a magnetic field is 
calculated in a relativistic Dirac-Hartree-Fock approach. In a mag netic field of 5 x 10 18 G in the 
center of the star, the overall equation of state is softer in comparison to the field-free case resulting 
in the reduction of maximum mass of the neutron star. 
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The matter density in the core of a neutron star could 
exceed up to a few times the nuclear matter saturation 
density. At such high density, it is expected that the 
quark degrees of freedom would be manifested. In fact, 
quark matter composed of comparable proportions of up, 
down and strange quarks has been conjectured to be 
the true ground state of QCD at finite baryon density. 
Therefore, at such high baryon density, a transition from 
nuclear matter to a stable quark matter is a possibility. 
Several authors J3]-^| have studied the effect of this phase 
transition on neutron star properties. 

The presence of strong magnetic fields in neutron stars 
might have interesting astrophysical implications. Large 
magnetic field B m ~ 10 14 G has been estimated at the 
surface of neutron stars |7|. On the other hand, in the 
core, the field may have been amplified considerably due 
to flux conservation from the original weak field of the 
progenitor during its core collapse. In fact, field as large 
as ~ 10 18 G in the core is predicted || using scalar virial 
theorem which is based on Newtonian gravity. At such 
high matter density, the effect of general relativity is sig- 
nificant and this gives rise to a very strong gravitational 
force H on the star. Consequently, the value of B m is 
expected to be further increased above 10 18 G. Because 
of highly conducting core, such a high field is frozen in 
jl0| and may not manifest at the surface. The energy 
of a charged particle changes significantly in the quan- 
tum limit if the magnetic field is comparable to or above 
a critical value Bm [fllf] , and the quantum effects are 
most pronounced when the particle moves in the lowest 
Landau level. The interaction of charged particles with 
strongly quantizing fields has been shown to modify the 
gross properties of matter on the surface ||,[l2| as well as 
in the core of neutron stars fll|[l4| ■ 

Theoretical studies of Fock (exchange) term, relevant 
to the star surface/crust, in intense magnetic field have 
been carried out using the simple Thomas-Fermi-Dirac 
model [lq ]. In this Letter, we investigate the composi- 
tion and structure of neutron stars with quark-hadron 
phase transition under the influence of strong magnetic 
fields in the Dirac-Hartree-Fock (DHF) approach within 
a mean-field approximation. This method is rather gen- 



eral, so it should be of correspondingly broad interest. 

We describe the calculation of the Fock term in pres- 
ence of a magnetic field in a general formalism within the 
a-uj model |l4] , ^6| . In a uniform magnetic field B m along 
z-axis, the Lagrangian is given by 

C = 4, [i 7M D" -m-g a a- g^u"] i> + \{d^a) 2 

i (d u cj v - duU^f + im 2 (w M ) 2 , (1) 
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in the usual notation fl4|. The general solution for pro- 
tons is ip(r) cx exp(—ie^t + ip y y + ip z z)fp ytPz (x), where 
the 4-component spinor at zero temperature f p yPz (x) is 
of the same form as in Ref. |14|, but with the single- 
particle Hartree energy e H replaced by the correspond- 
ing Hartree-Fock energy e HF . The form of the spinor in 
a magnetic field (see Ref. ) restricts the evaluation of 
the Fock term to strong fields such that only the zeroth 
Landau level, v = 0, is populated. The position depen- 
dent part can then be decoupled so that it reduces to the 
form 
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The DHF equation of pro- 



can then easily be written as 



\ot z P; 



+ (3{m + U H + U F )} u( Pz ) = e^ F u( Pz ), (3) 



where the effective mass m* = m + Uif + Ug (p z ), and 
u(p z ) is the momentum dependent part of the spinor. 
For B m 0, the Hartree contribution U H is given in 
Ref. p4| . The Fock term U F in general is given by M] 
(3U F (p z ) = /3U F (p z ) + U F ( Pz ) + a ■ pU F ( Pz ), where p 



is the unit vector along B n 
terms are given by 
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Here q vz = q z (l + U^(p z )/q z ) and J a = 
gl exp (l 2 a /2qB m ) f* /2 d9 sec 9 erf (l a sec 9/ ^/2q~B^) with 
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e ^~ F ) + m ai and erf(x) denotes 



the error function; g a s correspond to meson coupling 
constants with a = (<j,lu). The exchange contribution 
for B m — is given in Ref. Jl7],[l^]. It is straightforward 
to extend this formalism to calculate the exchange in- 
teraction for electrons in the outer crust of a strongly 
magnetized neutron star (ljjl^l and also for a magnetized 
pair Fermi gas p0| at finite temperature and baryon den- 
sity in exotic stellar objects and cosmology. To explore 
the quark-hadron phase transition, we demonstrate here 
the extension of the DHF approach to the calculation of 
the one-gluon exchange term in the quark phase. 

The pure quark phase consisting of u, d and s quarks 
interacting through one-gluon exchange in local charge 
neutral and /3-equilibrium conditions is described by the 
bag model Q . The DHF equations of motion for quarks 
in strong magnetic field can readily be obtained from 
Eq. (3) by dropping the cr-meson term and replacing 
the w-meson coupling by the quark-gluon coupling, i.e. 
g u — > {g/2)X 1 ab , where a(b) corresponds to the color 
charge of the outgoing(incoming) quark and As are SU(3) 
generators. The Hartree term (U H ) vanishes because the 
color symmetric combination (~ trA) does not couple to 
the gluons. The interaction energy density is then due to 
exchange term only, and this to order g 2 for each flavor 
with B m ^ is 



in Ref. H|, but the first term corresponds to those 
for quarks with df = 3 and the second term is that 
for electrons. The total energy density of the pure 
quark phase is then £ q — £^ + J^f £j' I/ ~° + £ m + B, 
where B is the bag constant and the magnetic field 
energy density is £ m — Bf n (nb / n ) / (8ir) . Since the 
higher order contributions of the density dependent field 
Bmirib/no) (given by Eq. (5)) to number densities 
and chemical potentials of different species are found 
to be negligible, the magnetic energy density and mag- 
netic pressure have been treated perturbatively. The 
pressure follows from the relation P q — YfHf n f + 
fi e n e — £ q , where Hf denotes the quark chemical poten- 
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tial and n f = (d f qfB m /2-K 2 ) YZ^a 9u {nf - 
is the quark density. The electron density is n 
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(eB m /2n 2 ) 9u {hI ~ m l,u) > and He its chemi- 

cal potential. The charge neutrality condition, Q q = 
Yfqfn-f — n e — 0, and the /3-equilibrium conditions, 
Hd = Ha + He = Hs, can be solved self-consistently to- 
gether with the effective masses at a fixed baryon number 
density n q = (n u + + n s )/3 to obtain the equation of 
state (EOS) for the deconfined phase. For the ease of nu- 
merical computation we, however, add here the one-gluon 
exchange term perturbatively to energy density and pres- 
sure. 

To describe pure hadronic matter consisting of neu- 
trons (n), protons (p) and electrons (e), we employ the 
linear a-uj-p model of Ref. ^lj in the relativistic Hartree 
approach. Neglecting the Fock term in hadron phase is 
quite justified as the Hartree energy grows like p F while 
the exchange energy behaves as p\ E^] , so the Fock cor- 
rection to the EOS is expected to be small |l8|]22|l . The 
EOS for this phase is obtained by solving self-consistently 
the effective mass in conjunction with the charge neutral- 
ity and /3-equilibrium conditions, Q h = n p — n e = and 
Hn = Hp + He at a fixed baryon number density . Here 
rii and Hi denote the number density and chemical po- 
tential; the subscript i refers to n, p and e. The energy 
density £ h = £o+£ m (£q is the contribution to the energy 
density due to nucleons and electrons ||l4|| ) and pressure 
P h in this phase are related by P h = UiHi — £ h ■ 

The mixed phase of hadrons and quarks comprising 
of two conserved charges, baryon number and electric 
charge is described following Glendenning The con- 
ditions of global charge neutrality and baryon number 
conservation are imposed through the relations xQ h + 
(1 — x)Q q — and n b — x n \ + (1 — x) n t, where x 
represents the fractional volume occupied by the hadron 
phase. Furthermore, the mixed phase satisfies the Gibbs' 
phase rules: Hp — + Hd and P h = P q . The total en- 
ergy density is £ = x£ h + (1 — X)£ q - 

In the present calculation the values of the dimension- 
less coupling constants for a, lo and p mesons determined 
by reproducing the nuclear matter properties at a satu- 
/ or e. The notation in Eq. (4) is same as that ration density of no = 0.16 fm -3 are adopted from Ref. 
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where qj, m*j and p F are the charge, effective mass and 
Fermi momentum of quark of fth flavor. The QCD cou- 
pling constant is defined by a c — g 2 /4tt. The general 
expression (for all v) for the total kinetic energy of the 
quark phase in a magnetic field is 
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TABLE I. The phase boundaries, u\ and u 2 , and central densities u c of stars with maximum masses M max /M0 with 
and without magnetic fields that undergo a quark-hadron phase transition with interacting quark phase. The corresponding 
quantities with non-interacting quark phase are shown in parentheses. The variation of magnetic field with density n& is given 
by Eq. (5) with f3 = 0.01, 7 = 3 for B = 5 x 10 18 G and B = 10 19 G. Calculations are performed for the hadronic Lagrangian 
of Zimanyi and Moszkowski (ZM) [21] and Serot and Walecka (SW) [16]. The bag constant is B = 250 MeV fm~ 3 and the 
nuclear matter saturation density no is 0.16(0.1484) fm -3 for ZM(SW) model. 



Hadronic Model 


B m (Gauss) 


ui = m/n 


u 2 = n 2 /n 


u c = n c /n 


M max /M 


ZM 




10 14 - 5 x 10 18 
10 14 - 10 19 


4.107(3.329) 
4.108(3.326) 
4.130(3.312) 


25.941(17.681) 
25.754(17.530) 
25.383(17.085) 


7.152(7.754) 
6.902(7.104) 
5.214(5.016) 


1.707(1.610) 
1.549(1.466) 
1.331(1.268) 


SW 




10 14 - 5 x 10 18 


2.158(1.953) 
2.159(1.945) 


9.523(7.837) 
9.374(7.769) 


4.259(4.012) 
3.252(3.554) 


2.594(2.279) 
2.342(2.017) 



H . The current masses of u and d quarks are taken as 
m u = m d = 5 MeV and m s = 150 MeV, and the QCD 
coupling constant is a c — 0.2. We consider the bag con- 
stant B = 250 MeV fm~ 3 which corresponds to the lower 
limit dictated by the requirement that, at low density, 
hadronic matter is the preferred phase. The variation of 
the magnetic field B m with density rib from the center to 
the surface of a star is parametrized by the form 

B m {n b /n ) = £C rf + B [1 - exp {-/3(n fc /n ) 7 }] , (5) 

where the parameters are chosen to be f3 = 0.01 and 
7 = 3. The maximum field prevailing at the center 
is taken as Bq = 5 x 10 18 G and the surface field is 
B^ ~ 10 14 G. The number of Landau levels populated 
for a given species is determined by the field B m and 
baryon density fll4| . Over the entire density range from 
the surface to the core, we find that the nucleonic energy 
density (pressure) dominates substantially over the mag- 
netic field energy density (pressure). Hence matter can 
sustain such a high field without injecting any instability 
in the corresponding EOS. 

With this parameter set, we show in Table I, the mixed 
phase boundaries of neutron star for B m = and for 
B m =/= with interacting quark phase (IQP). To examine 
the measure of the importance of the exchange contri- 
bution, the corresponding results for matter with non- 
interacting quark phase (NQP) are given in parenthe- 
ses. The onset of transition is at density n% = uiUq and 
the pure quark phase begins at density 712 = 1*2 1-0 • F° r 
matter with NQP and B m = 0, the boundaries are at 
u\ = 3.329 and U2 — 17.681. With the inclusion of inter- 
action, the EOS for the quark phase becomes softer and 
transition to a mixed phase is delayed to u\ = 4.107. At 
high density, as expected, the EOS for the quark phase 
is more softer resulting in a larger shift in 112 to 25.941 
compared to that in u\. As a consequence, the extent of 
the mixed phase is increased. With further inclusion of 
magnetic field Bo = 5 x 10 18 G, the situation is reversed. 
The EOS for the hadronic sector for B m ^ is found to 
be softer at low density and it turns out to be stiffer at 
high density compared to the field- free case O] . On the 



other hand, the EOS for the IQP with B m ^ is found 
to be considerably stiffer over the entire density range 
compared to those for B m = 0, with or without the ex- 
change term. Therefore, for Bq = 5 x 10 18 G, the onset 
of mixed phase occurs at 4.108. The stiff EOS for both 
hadronic and interacting quark sectors at Bo = 5 x 10 18 G 
reduce the upper boundary to a density of 1*2 = 25.754 
and thereby also the extent of the mixed phase. Since 
the EOS of NQP with B = 5 x 10 18 G is most stiff of all 
the cases considered above, the boundaries of the mixed 
phase are maximally shifted to the lowest baryon densi- 
ties at u\ — 3.326 and 1*2 = 17.530, and thus has the 
smallest mixed phase extent. 

For the field-free case, with the appearance of quarks, 
the neutron and electron abundances are found to de- 
crease since quark matter furnishes both baryon number 
and negative charge. With Bq = 5 x 10 18 G, apart from 
reducing the mixed phase extent, the magnetic field en- 
hances the electron fraction in the hadronic sector Jl4| . 
Because of charge neutrality condition n p = n e in the 
hadronic phase the proton fraction is also increased. The 
enhanced electron abundance persists even in the mixed 
phase. In the mixed phase, with B = 5 x 10 18 G, the 
u-quark abundance is found to be enhanced, while those 
of d and s quarks remains practically unaltered. 

The maximum masses of the stars M max , obtained by 
solving the TOV equation || are given in Table I for dif- 
ferent cases studied. With only nucleons and electrons, 
without any quark, the maximum masses of the stars are 
1.778M and 1.643M for B m = and B = 5 x 10 18 G. 
This is a manifestation of the softening of the EOS in 
magnetic field for these stars |14|| . However, the intro- 
duction of quarks, which softens the overall EOS, causes 
the maximum mass to be smaller. For B m = 0, the inclu- 
sion of one-gluon exchange in the quark sector causes a 
delayed appearance of the pure quark phase, and this re- 
sults in a stiffer overall EOS with a larger maximum mass 
of 1.707-Mq in contrast to matter with NQP. Due to fur- 
ther softening of the EOS by B = 5 x 10 18 G for matter 
with IQP, the maximum mass is reduced to a value of 
1.549M Q . The smallest M max = 1.466M is obtained 
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for the softest EOS with this value of Bo for matter with 
NQP. The central densities u c = n c /no of the maximum 
mass stars for all cases (see Table I) are less than u 2 and 
fall within the mixed quark-hadron phase. Therefore the 
presence of a pure quark phase is precluded. 

We have found that when the QCD coupling constant 
a c > 0.25, "strange matter" becomes negatively charged 
and requires positrons for overall charge neutrality and 
therefore would have catastrophic consequences ||]. 

Calculations are repeated with a higher value of Bq = 
10 19 G (see Table I) p3fl . Here the maximum masses of 
the stars for both IQP and NQP are found to be smaller 
than the observational lower limit of 1.44M imposed by 
the larger mass of the binary pulsar PSR 1913+16 
Therefore, the quark-hadron phase transition with the 
present hadronic EOS |2l[| sets a limit on the maximum 
value of the field that can be used for a stable system. 
To estimate the uncertainties in the parametrization of 
B m in Eq. (5), calculations are performed with IQP mat- 
ter for a rapidly increasing value of B m ((3 = 0.02 and 
7 = 3) and for a slowly increasing value of B m ((3 = 0.005 
and 7 = 3) taking Bq = 5 x 10 18 G. The respective max- 
imum masses are 1.464M and 1.664M Q . Apart from 
the uncertainties stemming from the parameters of the 
quark phase (B and a c ), the inadequate knowledge of 
the hadronic EOS at high densities relevant to a neutron 
star has important bearing on the phase boundaries in 
quark-hadron phase transition To explore this ef- 
fect, we have also performed calculations with a stiffcr 
hadronic EOS The phase boundaries for this EOS 

(see Table I) are found to be shifted at much lower densi- 
ties; the corresponding maximum masses are much larger 
than the observed values. Recently similar conclusions 
have been drawn j(| about the phase transition densi- 
ties in an effective field theoretic model including non 
linear scalar and vector meson interactions which soften 
the hadronic EOS. 

Neutron star cooling by URCA process may provide 
important informations about the interior constitution 
of the star. It is therefore of interest to see whether 
cooling by direct URCA process involving quarks can oc- 
cur for star with quark matter. The decay of d and s 
quarks are kinematically allowed (25) if they satisfy the 
respective inequality conditions pp—pp<pp< pp + p e F 
and p F — p e F < p s F < p F + pp. With B m — 0, direct 
URCA process involving u and d quarks is found to oc- 
cur for rib ~ no, while the latter relation is satisfied for 
rib ~ 5.3no. On the other hand, for Bo — 5 x 10 18 G 
(j3 = 0.01 and 7 = 3) the first inequality condition is sat- 
isfied for rif, ~ no, while the latter occurs for ~ 2.4no- 
Thus the threshold density for s quark decay depends 
sensitively on B m . Large neutrino luminosity and there- 
fore rapid cooling may serve as an observational means 
for the presence of strong magnetic field in the inner core 
of the star, and work in this direction is in progress. 



[3. 



[1] E. Witten, Phys. Rev. D30, 272 (1984). 
[2] E. Farhi and R.L. Jaffe, Phys. Rev. D30, 2379 (1984). 
J.C. Collins and M.J. Perry, Phys. Rev. Lett. 34, 1353 

(1975) ; G. Baym and S.A. Chin, Phys. Lett. B62, 241 

(1976) ; B. Freedman and L. McLerran, Phys. Rev. D17, 
1109 (1978); Strange Quark Matter in Physics and As- 
trophysics, Ed. J. Madsen and P. Haensel, Nucl. Phys. 
(Proc. Suppl.) B24 (1991); A. Drago, U. Tambini, and 
M. Hjorth- Jensen, Phys. Lett. B380, 13 (1996). 
N.K. Glendenning, Phys. Rev. D46, 1274 (1992). 
M. Prakash, J.R. Cooke, and J.M. Lattimer, Phys. Rev. 
D52, 661 (1995). 

H. Miiller and B.D. Serot, Nucl. Phys. A606, 508 (1996). 

C. Woltjer, Astrophys. J. 140, 1309 (1964); T.A. Mihara 
et al, Nature 346, 250 (1990); G. Chanmugam, Annu. 
Rev. Astron. Astrophys. 30, 143 (1992). 

D. Lai and S.L. Shapiro, Astrophys. J. 383, 745 (1991). 
S.L. Shapiro and S.A. Teukolsky, Black Holes, White 
Dwarfs and Neutron Stars (Wiley, New York, 1983). 
A. Muslimov and D. Page, Astrophys. J. Lett. 440, L77 
(1995); Astrophys. J. 458, 347 (1996). 
V. Canuto and J. Ventura, Fundam. Cosmic Phys. 2, 203 

(1977) . 

D.G Yakovlev and A.D. Kaminker, The Equation of 
States in Astrophysics, Ed. G. Chabrier and E. Schatz- 
man (Cambridge University Press, Cambridge, 1994), p. 
214. 

A. S. Vshivtsev and D.V. Serebryakova, JETP 79, 17 
(1994). 

S. Chakrabarty, D. Bandyopadhyay, and Subrata Pal, 
Phys. Rev. Lett. 78, 2898 (1997). 

I. Fushiki, E.H. Gudmundsson, C.J. Pethick, and J. Yn- 
gvason, Ann. Phys. (NY.) 216, 29 (1992) and references 
therein. 

B. D. Serot, Phys. Lett. B86, 146 (1979); B.D. Serot and 
J.D. Walecka, Adv. Nucl. Phys. 16, 1 (1986). 
Y.C. Leung, Physics of Dense matter (World Scientific, 
Singapore, 1984). 

M. Jaminon, C. Mahaux, and P. Rochus, Nucl. Phys. 
A365, 371 (1981); C.J. Horowitz and B.D. Serot, Nucl. 
Phys. A399, 529 (1983). 

B. Banerjee, D.H. Constantinescu, and P. Rehak, Phys. 
Rev. D10, 2384 (1974); A.M. Abrahams and S.L. 
Shapiro, Astrophys. J. 374, 652 (1991). 
J. Daicic, N.E. Frankel and V. Kowalenko, Phys. Rev. 
Lett. 71, 1779 (1993); Phys Rep. 237, 63 (1994). 
J. Zimanyi and S.A. Moszkowski, Phys. Rev. C42, 1416 
(1990). 

S.A. Chin, Ann. Phys. (NY.) 108, 301 (1977). 
Even with this high magnetic field of Bo = 10 19 G, the 
EOS is still found to be stable. 

J.M. Weisberg and J.H. Taylor, Phys. Rev. Lett. 52, 1348 
(1984). 

N. Iwamoto, Ann. Phys. (NY.) 141, 1 (1982). 



[8 
[9 

[10 
[11 

[12: 

[13; 
[i4; 
[is; 

[is; 
[17; 
[is; 

[19 

[20 

[21 

[22 
[23 

[24 

[25 



4 



